Influence of higher-order deformations in the ^^S + ^^^Er fusion reaction 
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The shape of the measured barrier distribution for the ''''S + ^®®Er reaction is analysed using the 
coupled-channels description. The ^^*Er nucleus is a good candidate to test current fusion models 
description of deformation since it has a large quadrupole deformation with an insignificant hexade- 
capole deformation. Coupling to weaker channels, the 2^ state in "^''S, the 3]~ state in ^''^Er, and 
the pair neutron transfer channel, all were found to have little influence on the barrier distribution. 
I A successful description of the barrier distribution was only obtained after the hexacontatetrapole 

»' ■ deformation term in ^®*Er {fi^) was included in the coupling scheme. However, a positive value for 

was needed where the macroscopic-microscopic model predicts a negative one. 



I. INTRODUCTION 



(N 

> 

' The influence of the internal structure of nuclei on the probability for heavy-ion fusion is well established. Coupling 
, of the relative motion of projectile and target nuclei to internal nuclear degrees of freedom produces a distribution of 
' fusion barriers that replaces the single fusion barrier which results when these couplings are not included The 
] degrees of freedom that affect fusion include the deformation of one or both reactants, the excitation of collective 

• surface vibrational modes, and single or multiplc-nucleon transfer channels. 

\ In general, collective modes involving large numbers of nucleons play the major role in determining the shape of 
^ ■ the fusion barrier distribution. Evidence for this can be seen in the measured barrier distribution for the system 
1) . 16q _|_ i44gj^ [HH], for example, where the lowest energy octupole state in ^^"^Sm accounted for a large part of the 
'"^ , structure in the experimental barrier distribution, whilst the transfer channels had only a relatively minor effect at 
^ ' energies in the barrier region. This difference is due to the weakness of the coupling of one- or two-particle transfer 
^ , channels relative to the coupling strength of the collective modes This does not imply, however, that there are no 

• significant effects of transfer on fusion. The effects of positive (effective) Q-M&hxe transfers, when present, are most 
\ evident in the enhancement of the fusion cross section at energies well below the average fusion barrier Bq, since they 
■ determine the energy of the lowest barrier |3[||,P . 

5— I ' Two examples of the sensitivity of the fusion barrier distribution to the effects of deformation are the ^^0-1- 
^ iS'^Sm and ^^O -I- ^^^W reactions Calculations which could reproduce the measured barrier distributions 

required not only a quadrupole deformation term but also inclusion of a hexadecapole term, where the sign of the 
latter was responsible for the difference in shapes of the two barrier distributions [compare the solid lines in Fig. 
|l|(a) and (b)]. Although the quadrupole and hexadecapole degrees of freedom are responsible for a large part of the 
barrier distribution's shape, there are small differences between the shape of the theoretical and experimental barrier 
distributions at centre-of-mass energies Ec.m. ^ Bq. It was shown in Ref. 1^ however, that if additional coupling 
to vibrational states and transfer channels were included in the theoretical model, then a better reproduction of the 
shape of the measured barrier distributions resulted. The calculations of Ref. |^ are shown by the dashed lines in Fig. 
|l(a) and (b) for the ^^''Sm and ^^^W targets, respectively. The deformation parameters obtained from these fits to 
the data were found to be in closer agreement with published (non-fusion) values [pT| than those calculations which 
excluded the additional couplings. 

However, as recognised in Refs. PJTo|, the addition of any weak couplings may result in an improvement in the 
reproduction of the shape of the measured barrier distribution. The resulting deformation parameters from the fit 
to the fusion data are in better agreement with non-fusion values since, to a large extent, the additional couplings 
in ^^^Sm and ^*^W mimic the changes in the deformation parameters required to reproduce the barrier distribution 
shape. 

In other measurements involving deformed target nuclei, a discrepancy between the theoretical and experimental 
barrier distributions in the region 0.94 < Ec.m./Bo < 1.0 MeV was found in possibly four different systems. Examples 
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are shown in Fig. |](c) and (d) for the ^^O + ^ssy ^^^^^ i2qi _|_ 232rp-|^ reactions, respectively. Is this discrepancy 
due to couphng to weaker channels or to some deficiency in the theoretical description of fusion involving deformed 
nuclei? Or are the measured fusion barrier distributions sensitive to higher-order deformations, such as the Pq 
deformation, presently absent from the analyses of fusion involving deformed nuclei? A recent paper [|l4| addressed 
this last question by examining the influence of the ffj deformation on fusion for the ^^O + ^^^Sm, ^^°W, and ^ssy 
systems. However, the disadvantage of this analysis [Q was that all the nuclei involved in the study have large /3j 
deformations, which meant it was difficult to isolate the effects due to the /3q term from the f3j contributions. 

In order to investigate the influence of higher-order deformations, such as the f3'^ term, a target nucleus with a pure 
prolate deformation should be chosen, to avoid distortion due to a hexadecapole component. The ^^^Er nucleus is a 
good candidate since it has a large prolate deformation and a hexadecapole deformation that is expected to be zero, 
or very small The /3j vibrational mode in ^^^Er is also weakly coupled compared to the ^ss^j target nucleus 

where this mode is significant. Additionally, it is desirable to choose a reaction that involves a projectile nucleus with 
a high charge because the width of the barrier distribution is proportional to Z1Z2, which produces a "magnification" 
of the coupling effects. 

The reaction ^''S -I- ^^*Er is suitable based on the above criteria. It was recently used |l^Jl^ in a measurement of 
the fusion cross sections, and the fusion barrier distribution, for input into transition state model calculations used to 
test flssion models at large angular momenta. In this paper, the fusion barrier distribution for ^"^S -I- ^^^Er from Refs. 
p7| , p^ is analysed to examine the role higher-order deformations play in heavy-ion fusion. Details of the experiment 
and results for the '^'*S -I- ^^^Er measurement have been published In Sec. || that follows, the sensitivity of the 
34g _|_ i68gj. Carrier distribution to the nuclear potential parameters is examined. Then the effects of deformation on 



the barrier distribution are calculated in Sec. Ill after checking the adequacy of the approximate treatment of the 
excitation energies of states in ^^^Er. In Sec. IV, the effects of coupling to the octupole vibration in ^^^Er, and states 
in ^''S is examined, followed by a discussion on the effects of transfer couplings, and the conclusion in Sec. 

II. NUCLEAR POTENTIAL PARAMETERS 

One uncertainty concerning the theoretical description of heavy-ion fusion is the true form of the nuclear potential 
between interacting nuclei in the absence of channel coupling. In fusion analyses, an energy-independent nuclear 
potential of Woods-Saxon form has often been used, with 

Vir) = \ , (1) 

1 + cxp[(r - roAp - roAji )/a] 

where Vq is the depth, rg is the radius parameter, and a is the diffuseness of the nuclear potential. One method for 
determining the parameters of the Woods-Saxon potential is the procedure described in Ref. where the fusion 
cross sections a have been fltted with a single barrier penetration model at energies « 10% above the average barrier 
where there is no longer any signiflcant barrier strength. The rationale behind this procedure is that at energies well 
above the average barrier the calculated cross sections are very sensitive to but relatively insensitive to the channel 
couplings, thus enabling an estimate of the "bare" nuclear potential. 

With the advent of precise fusion data, even this small sensitivity to the channel couplings at higher energies was 
measurable, and was dealt with usually in one of two ways. In calculations that included couplings with excitation 
energies smaller than the curvature of a single barrier fiujQ, a small increase in the fusion barrier was made in order to 
retain the quality of the fit to the high energy data obtained without the couplings. Where the excitation energies of 
states was greater than Tiloq, they were not included explicitly in the calculations, since these states only contribute 
to a renormalisation of the "bare" nuclear potential and have no effect on the shape of the barrier distribution . 

If these higher lying states are included in the calculations, care should be taken not be count their effects twice, or 
wrong conclusions about the position of average fusion barrier will be drawn. An example of this problem is shown in 
Ref. [ pO[ where the quoted "average barrier" is several MeV above the measured average fusion barrier as determined 
by their experimental fusion barrier distribution. 

The fit to the high-energy fusion cross-sections for '^^S -I- ^^*Er is shown by the solid line in Fig. H, where the 
quantity aEc.m. has been plotted on a linear scale for clarity. The parameters for the nuclear potential obtained from 
this fit are given in the second row of Table |, noting that the diffuseness obtained was a = 1.35 fm. The resulting 
average barrier is Bq = 123.1 MeV at a barrier radius of Rb = H-l fm. The diffuseness from this fit is much larger 
than the results from fits to elastic scattering data in Ref. where the fitted value for the diffuseness was found 
to be a = 0.63 fm. Large values for a, in the range a — 0.84-1.35 fm, have also been required to fit data in recent 
measurements |^,^,^ involving a range of projectile-target combinations. In Ref. |2^ an alternative form for the 
nuclear potential was suggested as a possible explanation for the large values of a. If the nuclear potential fell more 



2 



rapidly with increasing values of r than the rate suggested by the Woods-Saxon form of Eq. |l|, it may be possible to 
match the potential at large r required to fit the elastic scattering data, whilst retaining the fit to fusion cross-sections 
at the smaller values of r probed by fusion collisions. 

A clue as to the actual value for the diffuseness of the nuclear potential can be obtained from the slope of the fusion 
excitation function in the tunnelling regime. That is, at energies below the lowest barrier, where coupling effects no 
longer influence fusion. If the diffuseness obtained from fits to elastic scattering is appropriate for fusion, then the 
fusion excitation function will also have this slope provided the energy is sufficiently below the lowest barrier. If the 
fusion excitation function falls more rapidly with energy than a calculation with a = 0.63 fm, this means fusion takes 
place through a wider barrier (less penetrability), implying a larger value for a. In the present case, it is difficult to 
determine the slope of the fusion excitation function because of the large width of the '^'*S -I- ^^*Er barrier distribution, 
resulting from the target nucleus deformation, and the possible influence of positive Q value transfers. However, other 
experiments have recently being performed [ p3| , using reactions that have a suitably narrow barrier distribution, and 
preliminary analysis of this data does support a value for a significantly larger than the elastic scattering results. 

A calculation was performed to demonstrate the effect of using a diffuseness smaller than the fusion data require. A 
value of a = 0.65 fm was chosen, which overestimates the measured cross-sections at high energies (see the dashed line 
in Fig. 1^) . The nuclear potential parameters for this calculation are given in the third row of Table |. The potential 
parameters were chosen to keep the fusion barrier unchanged at Bq ~ 123.1 MeV, which meant that Rb was increased 
from 11.1 to 12.0 fm to compensate. The fusion barrier distribution resulting from a calculation with a — 0.65 fm 
and a quadrupole deformation of /3j — 0.338 |2j] only is shown by the dashed line in Fig. ^. The barrier distribution 
from this calculation also leads to a deterioration in the agreement with measured barrier distribution. 



III. EFFECTS OF DEFORMATION ONLY 



Having established that a smaller diffuseness parameter made the agreement between the calculated and experi- 
mental excitation function and barrier distribution worse, the effects of deformation are examined in this Section. 
Coupled-channels (CC) calculations were performed with the code CCDEGEN |2^, which is based on a version of 
the code CCFULL described in Ref. In CCDEGEN, the effects of deformation are calculated by coupHng to the 
ground-state rotational band of the deformed target nucleus, using the no-Coriolis approximation |^,^ with the 
excitation energies of the states in the ground-state band taken to be zero. The no-Coriolis approximation has been 



shown to be adequate for reactions involving heavy-ions pi 29 1. These approximations enable the CC equations to be 



decoupled, and the resulting eigenchannel equations, which correspond to fusion of the (inert) projectile and deformed 
target nucleus whose symmetry axis is orientated at an angle 6 with respect to the beam axis, are solved p^ ] to obtain 
the tunelling probability P,j{Ec.m., 0) for each partial wave Jh. The fusion cross section is then calculated using 

0'(£'c.m.) = TT A'^(2J+l)Pj(^;c.„0, (2) 
J 

where Pj (-Ec.m.) is the total tunnelling probability for each J averaged over all orientations, and is given by 

P.j{E,.n..)^\ Pj{E,,„,,,0) Sin 9d9. (3) 

Equation^ is exact for the classical situation, where the number of (degenerate) states in the rotor tends to infinity [^ . 
In actual rotational nuclei, the number of states is finite, and the integral in Eq. |^ is evaluated up to some maximum 
value of spin /max (where all spins are even) This is done using an n-point Gaussian integration formula, where 
n/2 = (/max + 2). For the calculations described in this work, for systems other than ■^^S + ^^^Er [see Fig. |l|(c) 
and (d)], a total of 5 states up to the 8+ state were included, which corresponds to a coupled-channels calculation 
with 5 channels. Including states with spins higher than the 8"*" state did not make an appreciable difference to the 
calculated fusion cross section. For the '^'^S -I- ^^^Er calculations, 6 states, up to the 10+ state, were needed to evaluate 
the integral to the desired accuracy. 

To test the zero excitation energy approximation used in these calculations, a comparison was made with a cal- 
culation which takes into account the finite excitation energies of the rotational states. The calculation shown by 
the dashed line in Fig. ^(a) was made with the CC code CCFULL which included coupling to all orders in the 
deformation parameter for the nuclear coupling potential, with the energy of the first 2+ in ^^^Er at 79.8 keV. The 
excitation energies of higher members of the rotational band were calculated according to the usual formula for a 
rotating rigid body. The solid line in Fig. |4|(a) is the barrier distribution from the CCDEGEN calculation, using the 
geometric description of Eq. ^ The parameters for the nuclear potential were identical in these two calculations (see 
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the second row of Table |). Both calculations were made with = 0.338 |2j] and included coupling to all orders in 
the deformation parameter for the nuclear coupling potential. The difference between the two calculations is small 
which can be attributed to the large deformation, or equivalently the low energy of the first excited state in ^^^Er, 
although the agreement is not perfect js^ . 

Having shown numerically that the zero excitation energy approximation in CCDEGEN is a good approximation to 
better than the accuracy of the data, it was used in all the fusion calculations that follow to test the effects of various 
deformation parameters. 

1. Quadrupole deformation 

The shape of the barrier distribution that results from consideration of the quadrupole deformation only (solid line 
in Fig. ^(a) with /3j = 0.338) cannot account fully for the shape of the measured barrier distribution. Although 
the area under the calculation (when a = 1.35 fm) matches to within 2-49{|^ the area under the measured barrier 
distribution, the latter peaks at a lower energy than the theoretical curve and has a more "triangular" shape. Varying 
the magnitude of does not improve the agreement. 

2. Hexadecapole deformation 

Although the hexadecapole deformation for ^^*Er is expected to be very small, a calculation was performed to 
ascertain the influence of this degree-of-freedom on the shape of the barrier distribution. The results of two CC 
calculations are shown by the broken lines in Fig. ^(b), one with /3j ~ 0.338 and f3j = +0.01 (dashed line) and the 
other with /3j = 0.338 and = —0.01 (dot-dashed hne). These two values for /3j span the range of hkely values 
for ^^^Er, as determined in Ref. where a value of /3j = —0.007 was estimated. Of the resulting range of shapes 
shown in Fig. ^(b), none lead to a significant improvement in the agreement with the data. 

3. Hexacontatetrapole deformation 

Since the effect of the hexadecapole deformation for the ^^^Er nucleus is small, what then is the effect of the 
hexacontatetrapole {(3'^) degree of freedom on fusion? Guidance as to the magnitude and sign of the /3j deformation 
comes from theory. The ground-state shapes for a large number of nuclei have been calculated by MoUer, Nix, 
Myers, and Swiatecki using a global macroscopic-microscopic model Q. They calculated a value of /3q = —0.025 
for ^^^Er |3^, which is a factor 3.5 times larger than the /?J estimate for the same nucleus. A CCDEGEN calculation 
with Pq = —0.025, in addition to the quadrupole deformation (/3j = 0.338), is shown by the dot-dashed line in Fig. 
^(c). The Pj deformation does have a significant effect on the shape of the barrier distribution, leading to more 
peaked shapes at each end of the barrier distribution. However, the inclusion of the negative hexacontatetrapole term 
worsens the agreement with the data. 

Another calculation was performed but this time with the opposite sign for the Pj deformation. This calculation, 
with P2 = 0.338 and Pq = -1-0.025, is represented by the solid line in Fig. ^(c). The inclusion of the Pq term with a 
positive instead of negative sign now improves the agreement with the measured barrier distribution, although all its 
features are still not completely reproduced. 

In the recent work of Ref. |Q, Rumin et al. investigated the influence of the Pq deformation on fusion for the 
^^O induced reactions on ^^''Sm, ^^^W, and ^^^U. In that work the authors concluded that the fits to the measured 
barrier distributions were improved by inclusion of the Pj terms, but in doing so obtained unphysically small Pj 
values for these actinides, contrary to their known hexadecapole deformations. It would be reasonable to conclude 
that the presence of the Pq term compensates for the large reduction in the Pj values. A contrasting approach is 
taken in this work, where the "known" values of /3j and Pj for ^^*Er are used rather than allowed to vary as free 
parameters in a flt to the barrier distribution. The basis for this approach is that the /3j value has been determined 
experimentally [Q, and the theoretical basis for Pj being close to zero is well established, since Er is midpoint the 
region of Z where Pj changes from a positive value (Sm and Gd, for example) to a negative one (Yb and Hf) |]l5| , |l6| . 



^The range depends upon the limits taken for evaluation of the integral. 
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An advantage the present study has over the analysis of Rumin et al. ||TJ] is that the Pj for ^^®Er nucleus is 
approximately seven times smaller than the /3j for the ^^^Sm and ^ssy nuclei (see Table |l|), implying there is very 
little influence of the hexadecapole deformation on the shape of the barrier distribution for ■^'^S + ^^^Er. In addition, 
the estimated |^ magnitude for the value for ^^^Er is around 1.5 times the estimate for ^ssy and 5 times the 
for ^^^Sm, increasing the possibility that fusion barrier distribution might be sensitive to the effects of this deformation 
parameter. When the deformation parameters for each nucleus are placed into the multipole expansion of the nuclear 
part of the coupling interaction, the hexacontatetrapole term for ^^*Er is 3.6 times larger than its hexadecapole term, 
compared to ^sSjj qq^qq where the hexacontatetrapole term is 3.4 times smaller than its hexadecapole term. So if the 
16q _|_ 238^ barrier distribution is sensitive to the presence of the /3j in ^^^U, then it may not be unexpected that 
the barrier distribution for ■^^S + -'^^^Er is sensitive to the Pq deformation in -"^^^Er. 

Although using the magnitude for Pg from MoUer's estimate produced a good description of the measured barrier 
distribution, the sign had to be inverted in order to achieve this agreement. In Ref. [p^, a similar problem arose, 
where the sign for Pq was also found to be opposite to theoretical predictions, after fits to the experimental barrier 
distributions for ^^O + 238^ ^^^^ i6q _^ i54g]^ were performed. An attempt was made to address this sign problem by 
re-fitting the barrier distributions with, in addition to the Pj, and Pj terms, coupling to the octupole vibration 
of 238u g^j^j^ 15'^Sm. This calculation then matched the theoretically predicted sign for Pq but, as discussed above, 
unphysically small values for Pj 10 were obtained. In the case of ■^'*S + ^^^Er, it was not possible to explain its sign 
problem by coupling to the 3j~ state in -'^^^Er, because this state couples so weakly (see Sec. IV below). 

To check that the positive sign required is not due to truncation of higher-order terms in the Coulomb coupling 
equation, another calculation was performed which included the (/3j)'^ terms of the Coulomb interaction. It was found 
that these terms had a very minor effect on the calculated barrier distribution for ^^S -I- ^^^Er, corresponding to a 
change by a "line width" in the shape of the barrier distribution, and could not explain why a positive sign for Pq 
was needed. 

Whether or not the measured barrier distribution for "^"^S -I- ^^^Er can determine the sign of the Pq deformation in 
i68gj. j-gjnains an open question. The macroscopic-microscopic calculations of Moller et al. predict strong negative 
Pe deformation in the N « 100, Z « 60 region, however it is worth noting the comment in Ref. [|32| that the behaviour 
of hexacontatetrapole deformation across the chart of the nuclides is less regular than that of the lower order, even 
(P2, P4) multipole distortions. 

Since the size of the deformation parameters extracted from the measured barrier distributions can depend upon 
the presence or absence of weaker couplings [0 , the effects of these channels are investigated in the next Section. 



IV. COUPLING TO WEAKER CHANNELS 



In this Section, coupling to channels other than the rotational coupling is considered, starting with the octupole 
vibration in '^^^Er. The dashed hue in Fig. |(a) is a CCDEGEN calculation with = 0.338, Pj = -1-0.025, and 
with the first 3~ state in "'^^^Er. The energy for this state is 1431 keV and the coupling strength was taken to be 
/3j = 0.064, obtained from the measured B{E3) T value Q with a radius for ^^^Er taken as 1.06^^^^ fm The 
barrier distribution with the octupole coupling is only slightly different to the calculation without it [solid line in 
Fig. ||(a)] implying that this state does not play a major role in the coupling scheme. This is expected since the 
coupling strength for the octupole state is very weak. Similarly, coupling to higher lying states, for example 7- and 
/3-vibrations, should have very little influence on the barrier distribution, as shown to be the case in Ref. |Q for ^^O 
induced reactions on the deformed ^^^Sm, ^^^W, and ^ssy targets. The effects of coupling to more exotic modes of 
excitation p4|, such double- /3 vibrations (2-phonon collective excitations), have not been examined. 



4- Coupling to states in S 

Up to this point in the analysis, '^^S has been treated as inert in the calculations. However, low lying states in ''^S 
will contribute to the channel coupling, and the size of their influence should be calculated. The barrier distribution 
shown by the solid line in Fig. ^(b) is a calculation which includes /3j and Pg deformation in the target, plus coupling 
to the 2^ state in '^''S, whose energy is 2.127 MeV with a coupling strength of P2 = 0.25. Also shown by the dashed 
line in Fig. |^(b) is the effect of coupling to the 2-phonon excitation in the projectile (only the [2^ (g)2^] 2-phonon state 
was included). Both these calculations have been shifted up in energy by 1.0 MeV to account for the adiabaticity of the 
projectile coupling |^,^. Coupling to the first excited state in the projectile does not dramatically alter the shape 
of the barrier distribution from the calculation which includes target excitations only. The single-phonon excitation 
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in the target produces a second peak in the barrier distribution centered around 133 MeV, but the uncertainty on 
rf^(-Ec.m.o')/d-Ec.m. is too large to identify such a feature. 

Coupling to the first excited state in ^^S has a small effect because the strong coupling from i^^Er produces a very 
broad barrier distribution. This means that the overall change in the shape of the barrier distribution is small, even 
though projectile coupling by itself produces a significant feature. This can be seen in Fig. ^, where a calculation 
with projectile coupling only is compared to a calculation which considers the target deformation only. 

Two other calculations were performed, again with = 0.338 and f3j — +0.025 for ^^^Er. The first included 
coupling to the 3|f state in '^^S at an excitation energy of 4.62 MeV and (3^ = 0.388 |^^. Coupling to this state 
only shifted the barrier distribution in energy, leaving its shape largely unchanged, as expected for a state which has 
such a large excitation energy |pl| , ^ . In the other calculation, inclusion of the 2f state in ^''S and = 0.338 and 
= —0.025 for ^^^Er, the theoretically predicted sign for the hexacontatetrapole deformation, produced a barrier 
distribution that in no way resembled the measured shape (calculation not shown). No reasonable combination of 
coupling to states in the projectile and the predicted value for f3j in ^^^Er could be found to reproduce the desired 
shape of the measured barrier distribution. 



5. Coupling to transfer channels 



Other weaker couplings not yet taken into account include transfer channels. The effective Q values, that is after 
taking into account the change in Z1Z2 of the transfer products (where relevant), for various transfer channels are 



shown in Table [II. Most of these transfer channels have a negative effective Q value, and the effect of their inclusion 
in a CC calculation would be to produce barriers with energies greater than the average fusion barrier [ p6[ . However, 
in the presence of the strong collective couplings in ^^*Er, the influence of the negative effective Q value transfers on 
the shape of t he bar rier distribution is likely to be small [g| , as was seen with the inclusion of the inelastic excitations 



S in Sec. [V 4. With the CC code used in this analysis it was not possible to calculate correctly the effects 
of all the above transfer channels. However, a calculation was performed including the 2n pickup channel, since 
the macroscopic form factor assumed for this pair transfer is likely to be close to the results from a microscopic 



treatment ||37|. Coupling to the pair transfer was made using the form factor 

Ft.a„(i?, 0) - -atil + piY2o + PlY,o + f^I^^o)^^^^^, (4) 

where R is the radius of the deformed target nucleus, is the orientation angle of the deformed (axially symmetric) 
target nucleus, and Yxq are the spherical harmonics. In Eq. ^ at is the transfer strength parameter. 

The results are shown by the dashed line in Fig. Q which, in addition to the /3j and /3q deformation, with /3j = 0.338 
and = +0.025, includes coupling to the 2n transfer channel with Q = +2.7 MeV and at = 0.2 MeV, the latter based 
on previous estimates of the pair transfer coupling strength [0j3^]. The effect of the positive Q 2n pickup channel is 
to redistribute some barrier strength in the barrier distribution, although the overall shape of the barrier distribution 
is not changed significantly. For comparison, the dot-dashed line in Fig. ^ shows the effect on the barrier distribution 
when the transfer coupling strength is reduced to at = 0.1 MeV. An increase in the transfer strength parameter 
resulted in an increase in the width of the barrier distribution, producing a distribution wider than the experimental 
one. A calculation with 2n pickup channel and the theoretically predicted sign for /3j resulted in barrier distribution 
that was even more peaked than that shown by the dot-dashed line in Fig. ^(c). Neither of these calculations are 
shown. 



As found in Sec. IV 4, no reasonable combination of weak coupling, in this case the 2n transfer channel, and the 
theoretically predicted sign for (3^ (ie negative sign) could mimic the shape of the barrier distribution obtained for a 
positive value of (3q . It seems unlikely that coupling to known additional transfer channels could change the shape of 
the ■^^S + ^^*Er barrier distribution to the extent the in ^^^Er does and thus account for the barrier distribution 
shape without resorting to the inclusion of /3j in the CC calculations. 



V. SUMMARY AND CONCLUSION 



Detailed coupled-channels calculations have been performed in an attempt to describe the shape of the recently 
measured barrier distribution |l^,|l^ for the "^"'S + ^^*Er reaction. The ^^^Er nucleus was chosen to investigate the 
disagreement between theory and measurement at energies 1 or 2 MeV below the average barrier observed in a range 
of reactions on deformed target nuclei. Since ^^^Er has a very small hexadecapole deformation, and because the 3j~ 
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state in ^^^Er couples very weakly, it was expected that the presence of the quadrupole deformation alone would allow 
for a stringent test of the fusion model. 

The calculated shape of the barrier distribution from a coupled channels calculation which included the quadrupole 
deformation only, failed to match the experimental "^"^S + i^^Er barrier distribution. Agreement with the experiment 
was improved significantly when the hexacontatetrapole deformation was included, with /3j = +0.025. However this 
value for /3j has a sign opposite to that predicted in the macroscopic-microscopic model of MoUer et al. |32| . 

When the theoretically predicted value for 0^ was used, this resulted in a double-peaked barrier distribution in 
contrast with the data. No combination of coupling to low-lying states in the projectile and/or the 2n transfer 
channel, when using the theoretically predicted value for could reproduce the shape of the experimental barrier 
distribution. The best reproduction was obtained with a positive value for [3'^ . It was argued that if such higher-order 
deformations were to be visible in the fusion barrier distribution, then the '^^^Er nucleus is likely to be one of the best 
candidates to observe their presence, since it has a very small hexadecapole deformation. 

Within the framework of current fusion models, the fusion barrier distribution for "'^S + ^^^Er apparently defines 
the deformation of ^^*Er as positive, in contrast with theoretical predictions. Further theoretical investigations are 
warranted, which could also examine the possible influence on fusion of other higher-order deformations not considered 
in this analysis. 
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TABLE I. Parameters for the real nuclear potential for ^*S + ^''^Er [see Eq. (1)]. 

Vo (MeV) ro (fm) a (fm) 

392.5 0.80 1.35 

292.1 1.10 0.65 



TABLE II. Summary of the 02 , Pl , and /3j deformation parameters for nuclei as indicated. 



Nucleus f3j 


^ Pi 


_ Pi 




'""Er +0.338 [ 
'^''Sm +0.28 1 
238u +0.27. 


m -0.007 

1 +0.05 
3 +0.05 


P -0.025 
1 -0.005 

-0.015 


32] 
K] 
~1 



TABLE III. Summary of the effective Q values Qoffj for various transfer channels for the reaction S + Er. 



Transfer products channel Qcff 

3b g ^ ^"'Er n pickup -0.79 

33 g ^ 169 ^ Stripping -5.4 

35 _^ 167 P pickup -7.4 

33p _|_ i69rpjj^ p stripping +0.74 

36(2-1 166jjp ^ pickup -6.2 

32p _,_ i70rp^ J Stripping -2.8 

36g _^ 166j,j. 2n pickup +2.7 

32g _,_ iTOgj. 2n stripping -6.8 
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FIG. 1. Fusion barrier distributions for a range of deformed nuclei in the reactions (a) ^^O + ^^"Sm M, (b) + 
§0, (c) + 238u 0, and (d) + ^^^Th 
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FIG. 2. Linear plot of aEcm. for the measured fusion cross-sections (solid points) for ^'*S + ^®®Er compared with two 
single-barrier penetration model calculations, one with a = 1.35 fm (solid line), and the other with a = 0.65 fm (dashed line). 
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FIG. 3. The measured fusion barrier distribution for S + Er |_8|. Here the second derivative was evaluated with a step 
length A_Ec.m. ~ 3.33 MeV (solid circles) or A_Ec.m. ~ 6.66 MeV (solid squares). Also shown are two CC calculations which 
are identical except for the nuclear potential parameters, the solid line results from the calculation made with a = 1.35 fm, and 
the dashed line with a — 0.65 fm. See text. 
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FIG. 4. (a) Test of the zero excitation energy approximation used in the CC code CCDEGEN (solid line) compared with 
an exact calculation (dashed line). See text for details, (b) Barrier distributions showing the effects of including the /3j 
deformation in addition to The solid line is a CCDEGEN calculation with plj = 0, the dashed line with /3j = +0.01, and 
the dot-dashed with /3j = —0.01. (c) Effects of the hexacontatetrapole deformation on the barrier distribution. Shown are 

CCDEGEN calculations with /3j = 0.338 and Pj^O (dashed hne), Pj = 0.338 and (Jj = +0.025 (solid line), and /?J = 0.338 
and f3j = —0.025 (dot-dashed line). Note that f3j is zero in each of the calculations shown in panel (c). 
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FIG. 5. (a) The effect on the barrier distribution when the 3]" state in '^^^Er is taken into account (dashed line), with 
02 = 0.338 and 0q = +0.025. Also shown is the same calculation without the 3J" state (solid line), (b) Barrier distributions 
for two calculations with 02 and /3j deformation plus coupling to the 2^ state in ^''S (solid line) and, in addition, coupling to 
the 2-phonon state in ^''S (dashed line). 
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FIG. 6. Comparison of the barrier distribution that results from coupling to the first 2"*" state in S only (dashed line) and 
target deformation only (solid line). 
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FIG. 7. Barrier distributions with 0^ = 0.338 and /3e = +0.025 and the 2n pickup transfer channel with at = 0.2 MeV 
(dashed hne), crt = 0.1 MeV (dot-dashed Une), and without the transfer channel (solid line). Note that none of the barrier 
distributions in this figure have been shifted in energy. 
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